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Combined recent data from cosmic-ray detectors and gamma-ray detectors have produced some surprising 
insights regarding the sources of ultrahigh-energy cosmic rays (UHECRs), magnetic fields inside and outside the 
Milky Way, and the universal photon backgrounds. The energy-dependent composition of UHECRs implies a 
non-negligible contribution of sources located in the Milky Way, such as past gamma-ray bursts that took place 
in our Galaxy. Extended halos of distant sources seen in the Fermi data imply that intergalactic magnetic fields 
have average strengths of the order of a femtogauss. Such relatively low magnetic fields imply that the protons 
from distant blazars generate a detectable flux of secondary gamma rays in their interactions with the photon 
background. A comparison with the data shows an excellent agreement of the secondary photons with the spectra 
of distant blazars observed by atmospheric Cherenkov telescopes. 



1. Introduction 

Several surprising insights have emerged from 
the recent data on ultrahigh-energy cosmic rays 
and very high energy gamma rays. 

Pierre Auger Observatory (PAO) has reported 
a steady increase of the mean nuclear mass with 
energy between 2 EeV and 35 EeV [IE]. This 
unexpected result is difficult to reconcile with the 
usually assumed extragalactic origin of UHECR. 
Indeed, most candidate sources of UHECRs, such 
as active galactic nuclei ( AGN) , are characterized 
by a relatively high photon density in the region 
of acceleration, which implies that only a small 
fraction of nuclei can survive and come out of 
the source. In addition, UHECR nuclei are prone 
to disintegration in photon backgrounds, which 
limits their propagation distances. 

However, it was recently pointed out that, if the 
sources are located inside the Milky Way Galaxy, 
diffusion in the galactic magnetic fields could pro- 
duce the spectrum and composition pattern con- 
sistent with the PAO data [3]. This brings in 
sharper focus the possibility that past GRBs in 
Milky Way Galaxy, which are expected to occur 
on the time scales shorter than the diffusion times 
of heavy nuclei, can be responsible for a large frac- 
tion of UHECR between 1 and 30 EeV. Further- 
more, if the cosmic rays of 10 18 — 10 19 eV are 



nuclei produced in the Milky Way, the effects of 
diffusion in turbulent Galactic micro-Gauss mag- 
netic fields can can also explain the lack of galac- 
tocentric anisotropy [3]. 

The lack of absorption features in the gamma- 
ray spectra of distant blazars can be explained 
by the production of secondary photons by pro- 
tons emitted by these blazars |4|5j . This expla- 
nation holds as long as the intergalactic magnetic 
fields (IGMF) are not strong enough to deflect the 
protons significantly. The recent measurement of 
IGMF using the Fermi data shows these fields 
have femtogauss strengths t3 , which supports the 
aforementioned interpretation. 

2. GRBs and unusual supernovae in the 
Milky Way 

There is a growing evidence that long GRBs are 
caused by a relatively rare type(s) of supernovae, 
while the short GRBs probably result from the 
coalescence of neutron stars with neutron stars or 
black holes. Compact star mergers undoubtedly 
take place in the Milky Way, and therefore short 
GRBs should occur in our Galaxy. 

Although there is some correlation of long 
GRBs with star- forming metal-poor galaxies [7], 
many long GRBs are observed in high-metallicity 
galaxies as well 8 9 10 , and therefore one expects 
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that long GRBs should occur in the Milky Way. 
Less powerful hypernovae, too weak to produce a 
GRB, but can still accelerate UHECR [IT], with 
a substantial fraction of nuclei |12ll3j . 

If the observed cosmic rays originate from past 
explosions in our own Galaxy, PAO results have 
a straightforward explanation [3] . 

GRBs have been proposed as the sources of ex- 
tragalactic UHECR [14115113] , and they have also 
been considered as possible Galactic sources [TBI 
117118] . It is believed that GRBs happen in the 
Milky Way at the rate of one per t GRB ~ 10 4 — 10 5 
years |19I20I21I22| . Such events have been linked 
to the observations of positrons [23 24 25 26 27 . 

If local sources, such as past GRBs in the Milky 
Way, produce a small fraction of heavy nuclei, the 
observed fraction of UHE nuclei is greatly am- 
plified by diffusion. This is because the galactic 
magnetic fields are strong enough to trap and con- 
tain nuclei but not protons with energies above 
EcV. This observation leads to a simple explana- 
tion of the composition trend observed by PAO. 

3. Effects of diffusion 

As illustrated in Fig. 1, diffusion depends on 
rigidity, and, therefore, the observed composi- 
tion can be altered by diffusion |3|28j . Changes 
in composition due to a magnetic fields have 
been discussed in connection with the spectral 
"knee" [28], and also for a transient source of 
UHECR [29]. The "knee" in the spectrum oc- 
curs at lower energies than those relevant PAO, 
and at higher energies the cosmic rays effectively 
probe the spectrum of magnetic fields on greater 
spatial scales, of the order of 0.1 kpc [30] , 

One can use a simple model [3] to show how dif- 
fusion affects the observed spectrum of the species 
"i" with different rigidities. Let us suppose that 
all species are produced with the same spectrum 
n (src) _ n ( slc ) £j-7 a t the source located in the 
center of Milky Way and examine the observed 
spectra altered by the energy dependent diffusion 
and by the trapping in the Galactic fields. In 
diffusive approximation, the transport inside the 
Galaxy can be described by the equation: 

^-V(AVn») + ^(6 i n i ) = 




Figure 1. For each species, there is a critical 
energy Eq^ for which the Larmor radius Ri is 
equal to the magnetic coherence length l c . For 
E <C Eq^, the mean free path of the diffusing par- 
ticle is I ~ Iq, and D^E) = Z c /3. For E > E 0<i , 
the particle is deflected only by a small angle 
9 ~ lo/Ri, and, after k deflections, the mean de- 
flection angle squared is 2 ~ k(l /Ri) 2 . The 
corresponding diffusion coefficient is Di{E) oc 
(^-) 2 , for E » Eo ti . 



Qi(E,r,t) + J2 I P ik (E,E')n k {E')dE'. (1) 

k 

Here Di(E,r,t) — Di{E) is the diffusion coef- 
ficient, which we will assume to be constant in 
space and time. The energy losses and all the 
interactions that change the particle energies are 
given by bi(E) and the kernel in the collision inte- 
gral Pik{E,E'). For energies below GZK cutoff, 
one can neglect the energy losses on the diffusion 
time scales. 

The diffusion coefficient D(E) depends primar- 
ily on the structure of the magnetic fields in the 
Galaxy. Let us assume that the magnetic field 
structure is comprised of uniform randomly ori- 
ented domains of radius lo with a constant field 
B in each domain. The density of such domains 
is N ~ lg 3 . The Larmor radius depends on the 
particle energy E and its electric charge qi = eZf. 

Ri = l (Jr^j > where E o,i = E Z h (2) 

The spatial energy spectrum of random mag- 
netic fields inferred from observations suggests 
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that B ~ 3^G on the 0.3 kpc spatial scales, and 
that there is a significant change at I = 1/k ~ 
0.1 — 0.5 kpc [30 - This can be understood theo- 
retically because the turbulent energy is injected 
into the interstellar medium by supernova explo- 
sions on the scales of order 0.1 kpc. This energy 
is transferred to smaller scales by direct cascade, 
and to larger scales by inverse cascade of magnetic 
helicity. Singlc-cell-size models favor ^0.1 kpc 
scales as well [50] . 

As explained in the caption of Fig. 1, diffu- 
sion occurs in two different regimes depending on 
whether the Larmor radius is small or large in 
comparison with the correlation length. As a re- 
sult, the diffusion coefficient changes its behavior 
dramatically at E = Eqj: 

D <w = \ ;*"np-*o w 

Here the two parameters < <5i,2 < 0.5 are dif- 
ferent from zero if the magnetic domains are not 
of the same size. The exact values of these param- 
eters depend on the power spectrum of turbulent 
magnetic fields. 

The approximate solution of the transport 
equation in our simple model yields 

Since diffusion depends on rigidity, the composi- 
tion becomes energy dependent. Indeed, at criti- 
cal energy -Eb.i, which is different for each nucleus, 
the solution ([5]) changes from oc E^ 1 to oc _E~ 7 ~ 2 
because of the change in Di(E), as discussed in 
the caption of Fig. 1. Since the change occurs at a 
rigidity-dependent critical energy Eqj — eEoZi, 
the larger nuclei lag behind the lighter nuclei in 
terms of the critical energy and the change in 
slope. If protons dominate for E < Eq, their flux 
drops dramatically for E > Eq, and the heav- 
ier nuclei dominate the flux. The higher the 
higher is the energy at which the species experi- 
ences a drop in flux. 

One can also understand the change in compo- 
sition by considering the time of diffusion across 
the halo is ti ~ R 2 /Di. The longer the particle 




Figure 2. UHECR spectra according to the model 
of Ref . [3] , shown here for different values of mag- 
netic fields (to illustrate the robustness of the 
model). The magnetic field was assumed to be 
~ 10/iG, coherent over lo = 100 pc domains (c/. 
B ~ 4/iG, l = 200 pc in Ref. [3 ). The power 
and the iron fraction were adjusted to fit the PAO 
data [3"T] . 



remains in the halo, the higher is the probability 
of its detection. At higher energies, the magnetic 
field's ability to delay the passage of the particle 
diminishes, and the density of such particles drops 
precipitously for E > E oi . Since Ei is propor- 
tional to the electric charge, the drop in the flux 
occurs at different energies for different species. 

The model of Ref. [3] gives a qualitative de- 
scription of the data. To reproduce the data more 
accurately, it must be improved. First, one should 
use a more realistic source population model. Sec- 
ond, one should include the coherent component 
of the Galactic magnetic field. Third, one should 
not assume that UHECR comprise only two types 
of particles, and one should include a realistic dis- 
tribution of nuclei. Finally, one should include 
the extragalactic component of UHECR produced 
by distant sources, such as active galactic nuclei 
(AGN) and GRBs (outside the Milky Way). 

A recent realization that very high energy 
gamma rays observed by Cherenkov telescopes 
from distant blazars are likely to be secondary 
photons produced in cosmic ray interactions 
along the line of sight lends further support to 
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the assumption that cosmic rays are copiously 
produced in AGN jets |4l5j . For energies E > 
3 x 10 19 eV, the energy losses due to photodisin- 
tegration, pion production, pair production and 
interactions with interstellar medium become im- 
portant and must be included. The propagation 
distance in the Galaxy exceeds 10 Mpc, so that 
the Galactic component should exhibit an analog 
of GZK suppression in the spectrum. The extra- 
galactic propagation can also affect the composi- 
tion around 10 18 eV [32]. 

Galactocentric anisotropy for a source distri- 
bution that traces the stellar counts in the Milky 
Way is small 3 . Although the anisotropy in pro- 
tons is large at high energies, their contribution 
to the total flux is small, so the total anisotropy 
was found to be < 10%, consistent with the ob- 
servations. The latest GRBs do not introduce a 
large degree of anisotropy, as it would be in the 
case of UHE protons, but they can create "hot 
spots" and clusters of events. 

The model of Ref. 3] makes an interesting pre- 
diction for the highest-energy cosmic rays. Just 
as the protons of the highest energies escape from 
our Galaxy, they should escape from the host 
galaxies of remote sources, such as AGN. There- 
fore, UHECR with E>3x 10 19 eV should cor- 
relate with the extragalactic sources. Moreover, 
these UHECR should be protons, not heavy nu- 
clei, since the nuclei are trapped in the host galax- 
ies. If and when the data will allow one to de- 
termine composition on a case-by-case basis, one 
can separate E > 3 x 10 19 eV events into pro- 
tons and nuclei and observe that the protons cor- 
relate with the nearby AGN. This prediction is 
one of the non-trivial tests of our model: at the 
highest energies the proton fraction should ex- 
ist and should correlate with known astrophysical 
sources, such as AGN. The microgauss magnetic 
fields in Milky Way cause relatively small deflec- 
tions for the highest-energy protons. As for the 
intergalactic magnetic fields, the first detection 
and a measurement of such fields, using the data 
from Fermi Gamma-ray Space Telescope, points 
to relatively weak, femtogauss field strengths [6] , 
which should not affect the protons significantly 
on their trajectories outside the clusters of galax- 
ies. 



If local, Galactic GRBs are the sources of UHE- 
CRs, the energy output in cosmic rays should be 
of the order of 10 46 erg per GRB. This is a much 
lower value than what would be required of ex- 
tragalactic GRBs to produce the same observable 
flux. Indeed, in our model the local halo has a 
much higher density of UHECR than intergalac- 
tic space, and so the overall power per volume 
is much smaller. The much higher energy out- 
put required from extragalactic GRBs [14115113] 
in UHECR has been a long-standing problem. 
The same issue does not arise in our case be- 
cause it seems quite reasonable that a hypernova 
or some other unusual supernova explosion would 
generate 10 46 erg of UHECR with energies above 
10 EeV. 



4. Intergalactic magnetic fields 

A recent analysis of the stacked images of 170 
AGN observed by Fermi Gamma-Ray Space Tele- 
scope has made possible the first measurement 
of the intergalactic magnetic fields [6]. The com- 
posite image of multiple AGN would be pointlike, 
broadened only by the instrumental effects, in the 
absence of IGMF. However, the images obtained 
from Fermi data exhibit halos consistent with the 
deflections of electromagnetic cascades by IGMF. 
Two alternative analyses, based on the pre-launch 
calibration and on the actual image of the Crab 
pulsar, yield consistent and statistically signifi- 
cant evidence that the halos are real, and that 
they cannot result from an instrumental effects. 
Furthermore, the halo parameters depend on the 
redshifts of the sources in a way that is consistent 
with the effects of IGMF. 

The field strengths of IGMF inferred from this 
analysis are of the order of a femtogauss. This 
result may open a new window on the particle 
physics in the early universe which can be respon- 
sible for generating the cosmological magnetic 
fields in a phase transition, inflationary reheating, 
or another out-of-equilibrium phenomenon 33 34, 
155155] . 
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5. Protons and secondary photons from 
extragalactic sources 

As long as IGMF have strengths as low as few 
femtogauss [B], the protons from distant blazars 
travel essentially along straight lines and produce 
secondary photons in their interactions with pho- 
ton backgrounds |4|5j . For nearby sources, the 
primary gamma rays produced at the source dom- 
inate the ACT signals. However, for more distant 
sources, the primary photons are filtered out by 
their interactions with EBL and the secondary 
photons can make up most of the observed sig- 
nal [315] . 

The secondary photons provide an excellent fit 
to the spectral shapes of gamma-ray signals from 
distant blazars [415] . Such fits are essentially one- 
parameter fits that depend only on the overall 
power of the source emitted in cosmic rays. The 
spectral slope of protons and the level of EBL 
do not have a strong effect on the spectrum of 
secondary photons. However, for the same pho- 
ton flux, the neutrino flux varies depending on 
the maximal energy E max to which the protons 
are accelerated. Indeed, there are two compet- 
ing processes that generate secondary photons: 
PIebl -> P7T° -> P77 and VlCMB pe + e~. For 
smaller £ max , a larger fraction of photons come 
from the hadronic channel, which is accompanied 
by production of neutrinos via pjebl — >■ nir^ fol- 
lowed by the decays of charged pions and the neu- 
tron. Neutrino observations can help determine 
this parameter [5]. 

The sources of secondary photons exhibit an 
unusual scaling law, which can be useful for pop- 
ulation studies and for planning the future ACT. 
The flux of protons scales with distance d as 
1/d 2 , but the rate of their interactions in the (op- 
tically thin) photon background is proportional 
to d. Therefore, the flux of secondary photons 
scales as 1/d, not 1/d 2 . This scaling with dis- 
tance is the reason why the secondary photons, 
whose flux is proportional to become increas- 
ingly more important for distant sources, while 
the primary photon flux falls off more rapidly, as 
exp{— d/\}/d 2 , where A is the attenuation length 
due to absorption. 



1Q„GeV 

10]°GeV 

10 GeV 




Figure 3. Photon (low energy) and neutrino (high 
energy) spectra [5] expected from an AGN at 
z = 0.14 (such as 1ES0229+200), normalized 
to HESS data points (shown) [37], for i? max = 
10 8 GeV, 10 10 GeV, and 10 n GeV shown by the 
solid, dashed, and dash-dotted lines, respectively. 



6. Conclusions 

Based on the recent data, one can make several 
remarkable inferences about the ultrahigh-energy 
cosmic rays and magnetic fields inside the Milky 
Way and in the intergalactic space. First, the 
energy dependent composition, with heavier nu- 
clei at high energy, points to a non-negligible con- 
tribution from Galactic sources J3J- Diffusion in 
turbulent Galactic magnetic field traps the nuclei 
more efficiently than protons, leading to an in- 
crease in the nuclear fraction up to the energy 
at which iron escapes (~ 30 EeV). At higher 
energies, the extragalactic protons should dom- 
inate the flux of UHECR, and their arrival direc- 
tions should correlate with locations of the known 
sources. Second, the recent measurement of in- 
tergalactic magnetic fields shows that the have 
femtogauss strengths [5J, which means that high- 
energy protons generate secondary photons well 
aligned with the line of sight. This allows one to 
study the cosmic ray output of AGN using atmo- 
spheric Cherenkov telescope observations of dis- 
tant blazars. Finally, if and when neutrino tele- 
scopes, such as IceCube [35] detect point sources, 
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one can learn about the cosmic-ray sources and 
photon backgrounds by comparing the neutrino 
flux to the photon flux. Neutrino and gamma-ray 
observations can help distinguish the local Galac- 
tic sources from extragalactic sources of UHE 
nuclei [39140141] . These inferences open excit- 
ing new opportunities for multi-messenger pho- 
ton, charged-particle, and neutrino astronomy. 
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